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ABSTRACT 

Using the code AUTOSTRUCTURE, extensive calculations of inner-shell atomic data 
have been made for the chemical elements He, C, N, O, Ne, Na, Mg, Al, Si, S, Ar, 
Ca, Cr, Mn, Fe and Ni. The results are used to obtain up-dated opacities from the 
Opacity Project, OP. A number of other improvements on earlier work have also been 
included. 

Rosseland-mean opacities from OP are compared with those from OPAL. Differ- 
ences of 5 to 10% occur. OP gives the 'Z-bump', at log(T) ~ 5.2, to be shifted to 
slightly higher temperatures. The opacities from OP, as functions of temperature and 
density, are smoother than those from OPAL. 

The accuracy of the integrations used to obtain mean opacities can depend on 
the frequency-mesh used. Tests involving variation of the numbers of frequency points 
show that for typical chemical mixtures the OP integrations are numerically correct 
to within 0.1%. 

The accuracy of the interpolations used to obtain mean opacities for any required 
values of temperature and density depend on the temperature-density meshes used. 
Extensive tests show that, for all cases of practical interest, the OP interpolations give 
results correct to better than 1%. 

Prior to a number of recent investigations which have indicated a need for down- 
ward revisions in the solar abundances of oxygen and other elements, there was good 
agreement between properties of the sun deduced from helioseismology and from stellar 
evolution models calculated using OPAL opacities. The revisions destroy that agree- 
ment. In a recent paper Bahcall et al. argue that the agreement would be restored if 
opacities for the regions of the sun with 2 x 10 6 <T<5x 10 6 K (0.7 to OARq) were 
larger than those given by OPAL by about 10%. In the region concerned, the present 
results from OP do not differ from those of OPAL by more than 2.5%. 

Key words: atomic process - radiative transfer - stars: interiors. 



1 INTRODUCTION 

First results for opacities from the Opacity Project (OP) 
were published by Seaton et al. (1994, to be referred to as 
Paper I) and were in good general agreement with those from 
the OPAL project (Rogers and Iglesias, 1992) except for re- 
gions of high densities and temperatures. In a more recent 
paper, Badnell and Seaton (2003, Paper II) confirmed a sug- 
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gestion made by Iglesias and Rogers (1995) that the discrep- 
ancies in those regions were due to the omission of important 
inner-shell processes in the OP work. Further comparisons 
between OPAL results and those from OP with inclusion 
of all important inner-shell processes are given by Seaton 
and Badnell (2004, Paper III) for the case of the 6-element 
mix (H, He, C, O, S and Fe) of Iglesias and Rogers (1995). 
In Papers II and III, the required inner-shell atomic data 
were computed using the code AUTOSTRUCTURE (Badnell, 
1986, 1997). Similar computations of the required inner-shell 
atomic data for all cosmically-abundant elements (He, C, N, 
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O, Ne, Na, Mg, Al, Si, S, Ar, Ca, Cr, Mn, Fe and Ni) have 
now been made. Some further improvements and up-dates 
of the OP work are also considered in the present paper. 



2 OUTER-SHELL ATOMIC DATA 

Most of the original OP outer-shell atomic data were com- 
puted using sophisticated R-matrix methods (see The Opac- 
ity Project Team, 1995, 1997). Ions are specified by nuclear 
charge Z and number of 'target' electrons NE: the total num- 
ber of electrons in an ion is equal to (NE+1). For iron ions 
with NE= 13 to 18, the R-matrix data were supplemented by 
extensive data calculated by Lynas-Gray et al. (1995) using 
the code superstructure (Eissner et al. , 1974). Data for 
the less-abundant iron group elements, Cr, Mn and Ni, were 
obtained from iron data using interpolation techniques de- 
scribed in Paper I. 

All original outer-shell atomic data were computed in 
LS coupling and allowance for fine-structure was included 
using methods described in Paper I. 

In Paper III the experiment was made of replacing all of 
the original iron data for NE = 13 to 18 by data computed in 
intermediate coupling using AUTOSTRUCTURE. Use of those 
new data did not make any major changes, which provided 
a good check on the earlier work, but there were some im- 
provements due to the inclusion of intercombination lines 
and some improvements in photoionization cross-sections. 
The new intermediate-coupling iron data for iron (NE = 13 
to 18) are used in the present work. 

As in Paper I, we use data from Kurucz (1988) for the 
first few ionization stages of iron-group elements. 



for all allowed I', I". The contributions from higher- n in 
were obtained by matching onto the results of 10 . 
The L-shell processes that were included are 

2l q 3l' p n"l" +hu^> 2l q ~ 1 3l' p n"l" + e" (3) 

and 

2l q 3l' p n"l" +hv^> 2l q ~ 1 3l' p n"l"n"'l"' -> X(NE) + e _ , (4) 

where 3l' p stands for for 3s s 3p*3d" with p — s + t + u. 
Calculations were made for q — 1 — 7 for p — and q = 8 for 
p = — 2. We used n' , n", n'" = 3 to 6 for photoexcitation 
and = 3 to 5 for photoionization, for all allowed I', I", I'". 
The M-shell processes that were included are 

3/VZ' + hv -» Zl q ~ x n'l' + e" (5) 

and 

3/VZ' + hv -» ?,l q ~ 1 n'l'n"l" -> X(NE) + e". (6) 

Calculations were made for q — 1 and 2 and n' , n" — 3 to 6. 
Calculations for q = 3 (NE=13) for Fe were made previously, 
and reported-on in Paper II. 

The range of p and q used in equations Q - © are 
the same as were used in Paper II for Fe, while those used 
previously for C, O and S were more restrictive (but suf- 
ficient for the 6-element mix). The same range of p and q 
has now been used for all elements. Further details of the 
calculations with AUTOSTRUCTURE may be found in Paper 
II. 

The calculation using AUTOSTRUCTURE provided totals 
of 11523 624 spectrum lines and 815 989 photo-ionization 
cross-sections for use in the opacity work. 



3 INNER-SHELL ATOMIC DATA 

All inner-shell atomic data were calculated using AU- 
TOSTRUCTURE employing both Russell-Saunders coupling 
(LS) and intermediate coupling (IC). In Paper II it was re- 
ported that, for the inner-shell work, the use of LS data and 
allowance for fine-structure using methods described in Pa- 
per I, did not give opacities differing significantly from those 
obtained using the IC atomic data. The LS inner-shell data 
were used in the present opacity work. 

New inner-shell data have been calculated for the ele- 
ments N, Ne, Na, Mg, Al, Si, Ar, Ca, Cr, Mn and Ni for NE 
— 1 — 12. The present data were generated by working iso- 
electronically. Data for He, C, O, S and Fe were reported-on 
in Paper II and were generated by working iso-nuclearly. For 
C, O and S, this has led to the calculation of some additional 
inner-shell data (see below). 

The K-shell processes that were included are of the 
form: 



ls"2Z p n'/' +hv-> ls q - 1 2l p n'l' + e" 
for photoionization and 

ls q 2l p n'l' +hv^ Is"' 1 2l p n'l'n"l" -> X(NE) + c 



(1) 



(2) 



for photoexcitation (-autoionization), where 2l p stands for 
2s s 2p' with p — s + t. Calculations were made for q = 1 for 
p = and q — 2 for p = - 7. We used n', n" = 2 to 5, 



4 THE FREQUENCY MESH 

Let u — hv/(kT) and let cr(u) be the monochromatic opacity 
cross-section per atom. The Rosseland-mean cross-section is 
or where 



1 

cm 

and 



■F(u) du 



F(u) = [15/(47r 4 )]u 4 exp(-?t)/[l - exp(-u)] 2 



(7) 



(8) 



The Rosseland-mean opacity per unit mass is kr = (Jr/^i 
where is the mean atomic weight. 

The function F(u) has a maximum value for u — 
3.8300 . . . and is small for u small and u large. It is not prac- 
ticable to use a mesh in u such as to resolve all spectral-line 
profiles, the number of points would be much too large. We 
therefore seek a mesh such that trapezoidal-rule integrations 
give results of acceptable accuracy. In Paper I a mesh was 
used with 10 4 points in the range 0.001 ^ ii ^ 20. 

Equation J7J may be replaced by 



r{u) 



dv 



where 



v(u) = / F(u) du 
Jo 



(9) 



(10) 
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Table 1. Element abundances, A(n) = log(A r (n)), relative to 
A(H) = 12. 



Elem. 


BP04 


A04 


Elcm. 


BP04 


A04 


H 


12.00 


12.00 


Si 


7.56 


7.51 


He 


10.92 


10.925 


S 


7.20 


7.14 


C 


8.52 


8.41 


At 


6.40 


6.18 


N 


7.92 


7.80 


Ca 


6.35 


6.31 





8.83 


8.66 


Cr 


5.69 


5.62 


Ne 


8.08 


7.84 


Mn 


5.53 


5.46 


Na 


6.32 


6.17 


Fe 


7.50 


7.44 


Mg 


7.58 


7.55 


Ni 


6.25 


6.18 


Al 


6.49 


6.37 









and Umax = V (u — > oo). From numerical integration, « m ax = 
1.0553 .... With a given total number of integration points, 
NTDT, and constant intervals in u or v, use of in place of 
10 gives more points in regions giving large contributions to 
the integrals. It was found that the use of © gave a signif- 
icant improvement in the accuracy of calculated Rosseland- 
means. 



5 ELEMENT ABUNDANCES 

Results of Paper I were for a solar mix referred to as S92, 
based on work of Anders and Grevesse (1989) together with 
some later revisions. Abundances used by Bahcall and Pin- 
sonneault (2004) in constructing a solar model, which we 
will refer to as BP04, are not very different from those of 
S92 or of GN93 (Grevesse and Noels 1993) used in much of 
the OPAL work. However, recent work (Asplund et al. 2004, 
and references cited therein) has suggested a need for sub- 
stantial revisions in solar abundances for oxygen and other 
elements. Dr N. Grevesse has kindly provided us with an 
updated table of recommended abundances, which we will 
refer to as A04, taking account of recent work by Asplund 
and others. Table 1 gives, for the 17 elements considered in 
the present work, the abundances from BP04 and A04: BP04 
gives X = 0.7394, Z = 0.0170 while A04 gives X = 0.7403, 
Z = 0.0123. 

A detailed discussion of recent work on the solar chem- 
ical composition is given by Asplund, Grevesse and Sauval 
(2005). 

6 CONTRIBUTIONS OF INNER-SHELL 
TRANSITIONS TO ROSSELAND-MEANS 

We use the variable 

R = p/T'i (11) 

introduced in the OPAL work, where p is mass-density and 
T 6 = 10" 6 x T with T in K. The behaviour of R in some 
stellar models is shown on Figs. 1 and 2 of Paper I. 

Fig. 1 shows, for the S92 solar mix, log(ftR) against 
log(T) for log(i?) = -1,-2,-3,-4,-5 and -6, both with 
and without inclusion of inner-shell transitions. It is seen 
that the inner-shell transitions contribute very little to the 
means at the lower densities, say log(R) < —4. The reason is 
that, at a lower density, ions in a given ionization stage will 



have maximum abundance only at a lower temperature: and 
at the lower temperatures features in the monochromatic 
opacities due to inner-shell transitions will appear at larger 
values of the frequency- variable u corresponding to regions 
where the weighting function F(u) in equation J7J is small. 



7 COMPARISONS WITH OPAL 

OPAL data can be obtained from the OPAL web site 1 for 
any required composition. The basic data are given as tables 
of log(ftR,) as functions of log(T) and log(-R). Codes are given 
for interpolations to any required values of log(T) and log(p) 
and further interpolations in X and Z. 

The OP calculations of monochromatic opacities for 
each chemical element are made on a mesh of values 
of log(T) and log(A r c ) where N c is electron density (the 
meshes used will be discussed further in Section |5J. A code 
mixv.f reads the monochromatic opacities and calculates 
Rosseland-means on the OP (T,N e ) mesh. A further code, 
OPFIT.F, provides interpolations to any required values of 
T and p. It includes a facility to produce tables in OPAL 
format. Such tables will be used in comparing results from 
OPAL and OP. 

Fig. 2 shows log(«;R.) against log(T) for log(i?) = 
-1, -2, -3, -4 and -6, from both OP and OPAL. The over- 
all agreement is seen to be good. The feature at log(T) ~ 5.2, 
usually known as the 'Z-bump', is mainly due to transitions 
in iron ions with NE= 13 to 18. It is seen that, compared to 
OPAL, the OP feature is shifted to slightly higher temper- 
atures. 

Use of the logarithmic scale for log(ftR,), as in Fig. 2, 
does not allow one to see the finer details of the level of 
agreement between the two calculations. Percentage differ- 
ences, (OP-OPAL), are shown on Fig. 3 for log(i?) = -1.0 
to —2.5, and on Fig. 4 for log(7?) = —3.0 to —4.5. It is seen 
that the general level of agreement is in the region of 5 to 
10%. An excursion to larger differences at log(T) ~ 5.5, ap- 
proaching 30% at log(7?) = —4.5, is due to the differences 
in the high-temperature wing of the Z-bump, as shown on 
Fig. 2. 

The differences of Figs. 3 and 4 have a somewhat ragged 
appearance which results from some lack of smoothness in 
the values of kr, at the level of 1 or 2%. Fig. 5 shows 
<91og(KR)/(91og(T) at constant log(_R), calculated using first 
differences, for \og(R) = —2. The data from OP are seen to 
be smoother than those from OPAL. Use of the u-mesh in 
the OP work (see Section |1J gives an improvement in the 
smoothness. 



8 THE SOLAR RADIATIVE INTERIOR 

Helioseismology provides remarkably accurate values for the 
depth, Rcz, of the solar convection zone (CZ). With the 
earlier estimates of solar element abundances, such as S92 
or GN93 (see section 6) , values of Rcz and other data from 
helioseismology were found to be in good agreement with 
results from solar models computed using OPAL opacities. 
However, with the new abundances it is found in two recent 

1 www-phys. llnl. gov /Research / OPAL / 



4 N. R. Badnell et al. 



Table 2. Contributions to ftp, at the test-point, log(T) 
log(JVe) = 23.0. A04 abundances. 



6.35, 



Table 4. Ionization fractions for oxygen at the test-point (NE= 
— 1 for fully-ionized, NE= for hydrogenic). 



Skip 




NE 


OP 


OPAL 


None 


15.29 


-1 


0.415 


0.374 


H, He 


10.51 





0.471 


0.524 


O 


10.52 


1 


0.109 


0.098 


Fe group 


11.03 


2 


0.005 


0.005 






3 


0.001 


0.000 






4 


0.000 





Table 3. Occupation probabilities and level populations for hy- 
drogenic oxygen at the test-point. 







OP 


OPAL 


n 


W{n) 


P0P(n) 


W(n) 


P0P(n) 


1 


1.000 


0.361 


1.000 


0.325 


2 


0.987 


0.048 


0.993 


0.044 


3 


0.813 


0.048 


0.954 


0.050 


4 


0.151 


0.013 


0.830 


0.063 


5 


0.009 


0.001 


0.289 


0.031 


6 


0.000 


0.000 


0.062 


0.009 


7 






0.009 


0.002 


8 






0.000 


0.000 



Table 5. Rosscland-means at the test point. 



Mixture 



OP 



OPAL % difference 



A04 15.29 14.99 +2.0 

A04 less iron-group 11.03 10.63 +3.7 

Difference 4.26 4.38 -2.8 



papers to be necessary to increase opacities in the vicinity 
of Rcz'- by 19% according to Basu and Antia (2004); and by 
21% according to Bahcall et al. (2004a). In a further detailed 
study of the problem, Bahcall et al. (2004b) show that there 
are also discrepancies for the solar profiles of sound speed 
and density and of helium abundance. They argue that all 
such discrepancies would be removed if the opacities were 
larger than those from OPAL by about 10% in the region of 
2 x 10 6 K < T < 5 x 10 6 K (0.7 to 0.4R Q ). 

In the solar model of Bahcall and Pinsonneault (2004) , 
referred to as BP04, the base of the convection zone is at 
log(T) ~ 6.34 (r ~ O.715-R ). In the CZ the abundances 
are independent of depth and the values used for the model 
were similar to those of the S92 mix (i.e. not taking ac- 
count of recently proposed revisions). Below the convection 
zone, log(T) > 6.34 (r < 0.715-Rq), abundances depend on 
the nuclear reactions which have taken place and vary with 
depth. Fig. 6 shows the Rosseland-mean opacities adopted 
by BP04, using OPAL data. In the CZ the exact values of 
opacity are not of great importance for the calculation of a 
model, since in that region the temperature gradient is no 
longer controlled by the opacity, but comparisons of different 
opacity calculations is still of interest. 

Fig. 7 shows percentage differences, (OP-OPAL), be- 
tween opacities from OP and those used for BP04. In the 
region below Rcz the relative abundances of metals are as- 
sumed to be independent of depth but the mass-fractions 
X, Y and Z are variable. We use a code mixz.f to obtain 
monochromatic opacities for the specified mixture of metals, 
and a code mxz.f for variable X and Z (with Y = 1 — X — Z). 
At the lowest temperatures (say log(T) < 4.5) kr is varying 
very rapidly as a function of T (see Fig. 6) and the differ- 
ences (OP-OPAL) are as large as 10%. Throughout the rest 
of the CZ the differences are never larger than about 5%. 
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8.1 Examination of the Rosseland mean near the 
base of the convection zone 

In the region at and below the base of the convection zone 
the mean opacities from OP and OPAL do not differ by more 
than 2.5% (see Figure 7). In view of the current interest in 
the opacity in that region, we consider the results in some 
greater detail. 

In the BP04 model the base of the convection zone 
is at log(T) = 6.338, log(p) = -0.735. We consider the 
Rosseland-mean opacity at a nearby test-point on the OP 
(T,N e ) mesh, log(T) = 6.35 and log(iV e ) = 23.0 giving 
log(p) = -0.715 for both BP04 and A04. Figure 8 shows 
log(cr(u)) against u — hv/{kT) calculated using BP04 abun- 
dances. The Rosseland mean at the test-point is kr = 19.35 
cm 2 g _1 with BP04 and 15.29 with A04. 

8.1.1 Contributing elements 

The code mixv.f used to calculate Rosseland means has a 
"skip" facility: an element marked as "skip" is included in 
calculating the mass density for a given electron density but 
its contribution to cr(u) is omitted. Table 2 shows the result 
of skipping: (a) hydrogen and helium; (b) oxygen; and (c) the 
iron-group elements. It is seen that, for all three cases, sub- 
stantial contributions are skipped. We shall consider those 
three cases in greater detail. 

8.1.2 The contribution from hydrogen and helium 

OPAL and OP results for H and He are discussed in section 
7.1 of Paper III. The agreement is generally close but there is 
a region in the vicinity of log(T) ~ 6 where OPAL has larger 
populations in the ground-states of H and of He + leading to 
larger OPAL opacities. At our test-point, kr(OPAL) for an 
H/He mixture (X ~ 0.7, Y ~ 0.3) is about 5% larger than 
kr(OP), corresponding to a difference of 2 or 3% for the 
A04 mix. 

8.1.3 Oxygen 

Oxygen makes an important contribution to «r in the vicin- 
ity of Rcz because it has a high abundance and because 
about one half of the oxygen is in a hydrogenic state (see 
Table 4) and the oxygen Lyman lines and continuum come 
in a region where the weighting function F(u) is large (see 
Figure 8) . The atomic data for the hydrogenic stage should 
be quite accurate and those for other high ionization stages 
should also be of good accuracy since, for a highly ionized 
system, the electron-nuclear potentials will be large com- 
pared with electron-electron potentials. The main uncer- 
tainties are likely to be in the level populations, as discussed 
in Paper III. 

We do not have OPAL occupation probabilities W(n) 
for oxygen at the test-point, but we do have them for the 
rather similar case of carbon at log(T) = 6.0, log(p) = —2.0. 
From Table 1 of paper II we see that, compared with OP, 
OPAL has larger occupation probabilities for more highly 
excited states. From that table we obtain the approximate 
relation 

Wopal(ti) ~ W O p(0.716 x n) (12) 



(carbon atlog(T) = 6.0, log(p) = -2) 

where we interpolate in n as a continuous variable. We use 
(12) to obtain estimates of Wo pal (n) for hydrogenic oxy- 
gen. Table 3 gives populations W(n) and populations POP(ra) 
for hydrogenic oxygen at the test-point, from both OP and 
OPAL, and Table 4 gives the corresponding ionization frac- 
tions, f ion(NE) . 

The Rosseland-mean cross-sections for oxygen are or = 
4.011 x 10~ 4 atomic units (a 2 ,) from OP and 3.721 x 10~ 4 
with our estimates of OPAL occupation probabilities. The 
corresponding opacities for the A04 mix are 15.75 and 15.16 
cm 2 g _1 . The larger excited-state W(n) from OPAL affect 
the Rosseland means in two ways. Firstly, as noted in Paper 
II, they give a reduction the ground-state population and 
hence a reduction in the strengths of the Lyman absorption 
features. Secondly, they give increases in profiles of lines hav- 
ing more highly excited states and reductions in continuum 
absorption for dissolved lines (see Appendix A of Paper I). 
Both effects lead to a reduction in Rosseland means. 

Our value of kr from OP is larger than that obtained 
using our estimates of Wo pal for oxygen by 4%. 



8.I.4 Other light elements 

It is shown in Paper III that, in the vicinity of Rcz, carbon 
and sulphur have behaviours similar to that for oxygen, with 
c(OP) being a few per cent larger than cr(OPAL). It can be 
expected that the other light elements will have a similar 
behaviour. 



8.1.5 The iron group 

It is seen from Table 2 that elements of the iron-group make 
significant contributions to kr at our test-point. In the vicin- 
ity of that point the dominant ionisation stages for iron are 
oxygen-like and fluorine- like. The R-matrix atomic data for 
the iron-group used by OP should be of rather good ac- 
curacy. Table 5 gives values of kr at the test-point from 
OP and OPAL. The OPAL values are obtained using tables 
from the OPAL web site and the OPAL interpolation codes 
xztrin21.f. Results are given for: 

(i) A04 abundances; 

(ii) A04 abundances less contributions from the iron 
group (number fractions are set to zero for iron group el- 
ements, other number fractions are unchanged except for H, 
for which the number fraction is increased to ensure particle 
conservation) . 

It is seen that the OP contribution from the iron group is 
2.8 % smaller than the corresponding OPAL contribution. 



8.1.6 Conclusion 

We recall that Bahcall et al. (2004b) require kr, to be 
10% larger than kr(OPAL) for 2 x 10 6 < T < 5 x 10 6 K 
(0.7 to 0.4_Rq). In that range we find that kr(OP) never 
differs from kr(OPAL) by more than 2.5%. The closeness 
of the agreement is partly fortuitous: compared with OP, 
OPAL has larger contributions from H and He and from 
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the iron-group, but smaller contributions from oxygen. How- 
ever, an increase to 10% larger than OPAL in the region 
2 x 10 6 < T < 5 x 10 6 does not appear to be very plausible. 



9 CHECKS ON ACCURACIES OF 

INTERPOLATIONS AND OF FREQUENCY 
MESH 

Temperature and density indices I and J are defined by 
log(T) = 0.025 x I, log(7V c ) = 0.25 x J. (13) 
Four meshes are used: 

(i) fine mesh, 'f, AI = AJ = 1; 

(ii) medium mesh, 'm', A7 = AJ = 2; 

(iii) coarse mesh, 'c', AI — AJ — 4; 

(iv) very coarse mesh, 'C, AI — AJ = 8. 

The number of frequency points, equally spaced in the vari- 
able v (see Section gj, is taken to be equal to NTOT. Nor- 
mal OP production work is done using the 'm' mesh with 
even values of I and J, and NT0T= 10000. OPAL produc- 
tion work is done with 10000 equally spaced points in the 
variable u — hv/(kT), 0.002 ^ it ^ 20. Checks on the accu- 
racies of our integrations over frequency and interpolations 
in T and p were made for a solar mix. 

Firstly, check calculations were made using the 'C mesh 
with NT0T= 10000 and 20000. It was found that use of the 
finer frequency mesh (NT0T= 20000) never led to a change 
in the Rosseland-mean of more than 0.1%. 

Then, 'C calculations were made for odd values of I 
and J, so that none of the 'C mesh-points coincided with 
'm' mesh-points. Using values of Rosseland-means for the 
'm' points, interpolations to the 'C points were made using 
the code OPFIT.F (Seaton, 1993). For log(T) > 3.73 and 
log(7?) ^ —1.0 all differences are less than 0.6%. Larger 
differences occur only in two regions: small log(T), where Kb. 
is varying rapidly as a function of T (see Fig. 1) — the worst 
case is 1.6% at (log(T) = 3.525, log(R) = -7.2); and large 
log(ii), close to the upper boundary for OPFIT interpolations 
(worst cases are 3.1% at (log(T) = 6.525, log(i?) = -0.530) 
and 5.8% at (log(T) = 7.25, log(-R) = -0.531)). 

Bahcall et al. (2004) report that the use of OPAL opac- 
ities and different interpolation schemes yield opacity values 
in the solar-centre region differing by as much as 4%. In view 
of the importance of obtaining accurate opacities for that re- 
gion, we have made further checks. For 6.0 < log(T) ^ 8.0 
we have made 'm' mesh calculations with NTOT equal both 
10000 and 30000. For all densities considered in the OP work 
(log(N e ) = 15.5 to 23.5 for log(T) = 6.0, log(A" c ) = 22.0 
to 29.5 for log(T) = 8.0) we find the largest differences 
in single-element Rosseland-means, for NT0T= 10000 and 
30000, to be about 4%. Such comparatively large differences 
occur only for cases with deep minima in the monochro- 
matic opacities. For typical mixtures such minima get filled- 
in by contributions from other elements and the sensitiv- 
ities to NTOT become much smaller. For a solar-type mix, 
6.0 ^ log(T) ^ 8.0 and all densities, the differences between 
kr for NT0T= 10000 and 30000 are never larger than 0.1%. 

As a further check on interpolations we have made cal- 
culations using both the 'm' and 'P meshes with NT0T= 
10000. We used OPFIT.F to make calculations for a model 



with a fine mesh in log(T) and densities such as to give 
Iog(-R) = —1.75, roughly corresponding to a solar-centre 
model. For 6.2 ^ log(T) 7.5, values of kr from the 'm' 
and 'P calculation never differ by more than 0.4%: for the 
vicinity of the base of the convection zone they never differ 
by more than 0.1%. 



10 DATA ARCHIVES 

10.1 Outer-shell atomic data 

The OP outer-shell atomic data are archived in TOPbase, 
which is part of the OP archive at the CDS (Centre de 
Donnees de Strasbourg) 2 . The original version of TOPbase 
is described by Cunto et al. (1993). 

10.2 Atomic data from AUTOSTRUCTURE 

Both inner- and outer-shell data from autostructure have 
been archived according to the Atomic Data and Analy- 
sis Structure (ADAS) ad}38 and adf39 data formats (Sum- 
mers 2003). We have included LS and IC final-state re- 
solved photoexcitation-autoionization and photoionization 
data, summed over final channel angular momenta for both 
inner- and outer-shells. This will enable the data to be used 
in the future for other purposes such as the collisional- 
radiative modelling of photoionized plasmas. Further details 
of the data archival may be found in Paper II. 

10.3 Monochromatic opacities and computer 
codes 

Codes for the calculation of Rosseland-mean opacities and 
radiative accelerations for any chemical mixtures, tempera- 
tures and densities are described by Seaton (2005). A tar file 
containing the codes and all required monochomatic opac- 
ities is available for download from the TOPbase web site 
at CDS. Since the file size is a little under 700Mb, it is also 
available on a CD 3 . 

10.4 OPserver 

Rosseland mean opacities and their derivatives can also be 
computed online for any required chemical mixture using 
the OPserver at the Ohio Supercomputer Center, Columbus, 
Ohio, USA (Mendoza et al. , 2005). It can be accessed via 
the TOPbase web site at CDS. 



11 SUMMARY 

OP opacities have been up-dated by inclusion of all con- 
tributing inner-shell processes and some further refinements 
have been made. 

Comparisons with OPAL show agreement generally to 
within 5 to 10%. 

Using earlier estimates of metal abundances, good 
agreement was obtained between computed solar models 

2 http : //cdsweb . u-strasbg . f r/topbase 

3 Requests to: Claude. ZeippenSobspm.fr 
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and results obtained from helioseismology for Rcz and other 
data. With recent downward revisions in metal abundances 
that agreement is destroyed. It is shown by Bahcall et 
al. (2004b) that agreement would be restored if opacities 
in the region just below the base of the convection zone, 
2 x 10 6 < T < 5 x 10 6 K (0.7 to OARq), were 10% larger 
than those from OPAL . In that region we find kr(OP) never 
to be larger than kr(OPAL) by more than 2.5%. Results for 
opacities in the region near the base of the convection zone 
have been discussed in some detail. 

Studies have been made of the accuracies of the OP in- 
tegrations over frequency, used to obtain Rosseland means, 
and of the OP procedures used for interpolations in temper- 
ature and density. It is concluded that the OP Rosseland- 
means should be numerically correct to within 1% or better. 
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Figure 1. Rosseland-mean opacities from OP for S92 mix, 
with and without inner-shell contributions. 
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Figure 2. Rosseland-mean opacities from OP and OPAL 
for the S92 mix. 
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Figure 5. Values of 91og(K R )/91og(T) (at constant 7?) ap- 
proximated using first derivatives. From OP and OPAL for 
log(fl) = -2. 
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Figure 7. Percentage differences, (OP-OPAL), between 
opacities for a solar model from Bahcall and Pinsonneault 
(2004). The convection zone is the region with log(T) < 6.34. 
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10 

u=hv/(kT) 

Figure 8. The monochromatic opacity for the BP04 mix at 
log(T) = 6.35, log(iVe) = 23.0. Cross-section, a, in atomic 
units, tto§. 



